Spin-Dependent Transport of Electrons in a Shuttle Structure 
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We consider "shuttling" of spin-polarized electrons between two magnetic elec- 
trodes (half-metals) by a movable dot with a single electronic level. If the magneti- 
zation of the electrodes is antiparallel we show that the transmittance of the system 
can be changed by orders of magnitude if an external magnetic field, perpendicular 
to the polarization of the electronic spins, is applied. A giant magnetotransmittance 
effect can be achieved for weak external fields of order 1 H- 10 Oe. 

PACS numbers: 72.25.Hg, 73.43. Jn., 73.61.Ey, 72.50.Bb 

INTRODUCTION 

Metal-organic nanocomposite materials are interesting from the point of view of the 
"bottom- up" approach to building future electronic devices. The ability of the organic parts 
of the composite materials to identify and latch on to other organic molecules is the basis 
for the possible self assembly of nanoscale devices, while the metallic components provide 
mechanical robustness and improve the electrical conductance. 

Such composite materials are heteroelastic in the sense that the mechanical rigidity of 
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the organic and metallic components are very different. This allows for a special type of 
deformation, where hard metallic components embedded in a soft organic matrix can be 
rearranged in space at a low deformation energy cost associated with stretching and com- 
pressing the soft matrix. Strong Coulomb forces, due to accumulation of electronic charge 
in embedded nanoscale metallic particles, can be a source of such mechanical deformations. 
This leads to a scenario where the transport of electric charge, possibly due to tunnelling of 
electrons between metal particles, becomes a complex nano-electromechanical phenomenon, 
involving an interplay of electronic and mechanical degrees of freedom Such an interplay 
can lead to new physics, as was recently demonstrated theoretically for the simplest possible 
structure — a Nanoelectromechanical single-electron transistor. The electromechanical in- 
stability predicted to occur in this device at large enough bias voltage was shown to provide 
a new mechanism of charge transport Q]. This mechanism can be viewed as a "shuttling" 
of single electrons by a metallic island — a Coulomb dot — suspended between two metal 
electrodes. The predicted instabilily leads to a periodic motion of the island between the 
electrodes shuttling charge from one to the other. 

The shuttle instability appears to be a rather general phenomenon. It has, e.g., been 
shown to occur even for extremely small suspended metallic particles (or molecules) for 



which the coherent quantum dynamics o 



the^ tunnelling electrons P| or even the quantum 



dynamics of the mechanical vibration j^|5|,|6|,|7[ become essential. Nanomechanical transport 
of electronic charge can, however, occur without any such instability, e.g., in an externally 
driven device containing a cantilever vibrating at frequencies of order 100 MHz. A small 
metallic island attached to the tip of the vibratirig cantilever may shuttle electrons between 



metallic leads as has recently been demonstrated 



I . Further experiments with magnetic and 
superconducting externally driven shuttles as suggested in |9| , seem to be a natural extension 
of this work. Fullerene-based nanomechanical structures [lO| are also of considerable interest. 

The possibility to place transition-metal atoms or ions inside organic molecules intro- 
duces a new " magnetic" degree of freedom that allows the electronic spins to be coupled to 



mechanical and charge degrees of freedom 



12| |. By manipulating the interaction between 



the spin and external magnetic fields and/or the internal interaction in magnetic materials, 
spin-controlled nanoelectromechanics may be achieved. An inverse phenomenon — nanome- 
chanical manipulation of nanomagnets — was suggested earlier in A magnetic field, 
by inducing the spin of electrons to rotate (precess) at a certain frequency, provides a clock 
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for studying the shuttle dynamics and a basis for a dc spectroscopy of the corresponding 
nanomechanical vibrations. 

A particularly interesting situation arises when electrons are shuttled between electrodes 
that are half-metals. A half-metal is a material that not only has a net magnetization as 
do ferromagnets, but all the electrons are in the same spin state — the material is fully 
spin-polarized. Examples of such materials can be found among the perovskite maganese 
oxides, a class of materials that show an intrinsic, so called "colossal magnetoresistance" 
effect at high magnetic fields (of order 10-100 kOe) [isj ]. 

A large magnetoresistance effect at lower magnetic fields has been observed in layered 
tunnel structures where two thin perovskite manganese oxide films are separated by a tunnel 
barrier |l^, 14, Here the spin polarization of electronic states crucially affects the 

tunnelling between the magnetic electrodes. This is because electrons that can be extracted 
from the source electrode have there spins aligned in a definite direction, while electrons 
that can be injected into the drain electrode must also have there spins aligned — possibly 
in a different direction. Clearly the tunnelling probability and hence the resistance must be 
strongly dependent on the relative orientation of the magnetization of the two electrodes. 
An external magnetic field aligns the magnetization direction of the two films at different 
field strengths, so that the relative magnetization can be changed between high- and low 



resistance configurations. A change in the resistance of trilayer devices bv 



actors of order 



2-5 have in this way been induced by magnetic fields of order 200 Oe jl4, QQ. The 
required field strength is determined by the coercivities of the magnetic layers. This makes 
it difficult to use a tunneling device of the described type for sensing very low magnetic 
fields. In this paper we propose a new functional principle — spin-dependent shuttling of 
electrons — for low-magnetic field sensing purposes. We will show that this principle can 
lead to a giant magnetoresistance effect in external fields as low as 1-10 Oe. 

The new idea which we propose to pursue is to use the external magnetic field to manipu- 
late the spin of shuttled electrons rather than the magnetization of the leads. The possibility 
to "trap" electrons on a nanomechanical shuttle (decoupled from the magnetic leads) during 
quite a long time on the scale of the time it takes an electron to tunnel on/off the shuttle 
makes it possible for even a weak external field to rotate the electron's spin to a significant 
degree. Such a rotation allows the spin of an electron, loaded onto the shuttle from the spin- 
polarized source electrode, to be reoriented in order to allow the electron finally to tunnel 
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FIG. 1: Schematic view of the nanomechanical GMR device: a movable dot with a single electron 
level couples to the leads due to tunnelling of electrons, described by the tunnelling probability 
amplitudes TL^ji{t)), and due to the exchange interaction whose strength is denoted by JL,R{t). 
An external magnetic field H is oriented perpendicular to the direction of the magnetization in the 
leads (arrows). 

from the shuttle to the spin-polarized drain lead. As we will show below, the magnetic field 
induced spin-rotation of shuttled electrons is a very sensitive nanomechanical mechanism 
for a giant magnetoresistance (GMR) effect. 



FORMULATION OF THE PROBLEM. GENERAL EXPRESSION FOR THE 

CURRENT 

A schematic view of the nanomechanical GMR device to be considered is presented in 
Fig. 1. Two fully spin-polarized magnets with fully spin-polarized electrons serve as source 
and drain electrodes in a tunnelling device. In this paper we will consider the situation when 
the electrodes have exactly opposite polarization. A mechanically movable quantum dot 
(described by a time-dependent displacement x{t)), where a single energy level is available 
for electrons, performs forced harmonic oscillations with period T — 2t:/ ui between the leads. 
The external magnetic field is perpendicular to the orientation of the magnetization in both 
leads. 
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The Hamiltonian that governs the dynamical evolution of the system is 

n{t) = eo{a\a^ + a\ai) + ^(eaaj^^ia^L + £aai^Raa,R) (1) 

a 

— JL(t)(a|a| — ajflj) — jR{t){a^^ai — aja-f) — {g^H/2){a^.^ai + a|a|) 

a a 

where a}^ {aa^^R)) are the creation (annihilation) operators of electrons with the energy 
Sa on the left (right) lead (we have suppressed the spin indices for the electronic states in the 
leads due to the assumption of full spin polarization), 01^^(0-^(1)) are the creation (annihi- 
lation) operators on the dot, Eq is the energy of the on-dot level, JL{R)(t) = JL{R){xit)) 
are the exchange interactions between the on-grain electron and the left (right) lead, 
^L{R)(t) = ^L{R){x(t)) are the tunnel coupling amplitudes, g is the gyromagnetic ratio 
and /i is the Bohr magneton. 

The single-electron density matrix describing electronic transport between the leads may 
be presented in the form: 

a a 

Here |^"'-^) are single-electron states that obey the time-dependent Shrodinger {h = 1) 
equation with a Hamiltonian given by Eq. (^. The initial condition has the form 

^ _oo)) = \a,L{R))exp{-ieJ) , 

where \a,L{R)) is a single-electron state on the left (right) lead with energy 

We will suppose that the internal relaxation in the leads is fast enough to lead to equi- 
librium distributions of the electrons. This means that Wa,L{R) = fi^a T ^/2) (where f{e) 
is the Fermi distribution function) and V is applied voltage. 

The problem at hand is greatly simplified if one considers the large bias-voltage limit 

I V^-eo |>^^ALx, (3) 

where ly is the density of states on the leads. The restriction (jH} does not allow us to consider 
a narrow transition region of voltages from the zero-current regime at < ~ ^^^max 
the fully transmissive (in the absence of spin polarization effects) regime at V > eo + uA"^^^. 
However, it covers the in practise most important case when the fully transmissive junction 
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is strongly affected by electronic spin-polarization. Therefore, in our further considerations 
we will take Wa,L — 1, 'Wa,R — and Eq — 0. 

We will calculate the average current, /, through the system from the relation 



I 

ONr 
dt 



= dtTr{pj}, (4) 



ie[H, Nr] = ieTR{t) ^{a\aa,R - al j^ai) , 



where N is the electron number operator for the right lead, N — R'^a,R- 
In general, the state |\E'°'^) can be expressed as 

=c^(t)| T) + crWI i) + Y.^c-£^t)\(5,L) + cf{t)\(3,R)), (5) 

Thus the problem is reduced to determining the coefficients c^'^^^ and c'^s^-^y 
At this point it is convenient to introduce the bi-vectors 

so that the coefficients c^^^^ can be expressed as (see Appendix 1) 

c'^/ ^-i f (iiV^/'(*-*')T«(t')(e2,c"(0)- 

Here (a, b) is the inner product of two bi-vectors. As shown in Appendix 1, by using the 
wide band approximation {i.e. by taking the electron density of states in the leads v to be 
constant) the equation for the bi-vectors c° takes the form 

i^-^^mc^+nt). (6) 

Here P(i) = TL{t)e-''''^ei and the matrix R{t) is 

\ -gfxH/2 J{t) - irR{t)/2 J 
where J{t) = J hit) — JRif) and TLi^R-^it) = 27rz/A^^^(t) is the level width. 
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The formal solution of Eq. (joj) can be written in the form 

c"(t) = [ dt'L{t,t')r{t'), (8) 

J — oo 

where the "evolution" operator L{t, t'), {L{t, t) = /), is defined as the solution of the equation 

^^^%^ = mLM, (9) 

and obeys the multiplicative and periodicity properties, 

L{t,t') = L{t,t")L{t",t'), L{t + T,t' + T) = L{t,t'). (10) 
Using Eq. (jHl) together with Eq. (0}, one can write the average current on the form 

J = ^ rdtVnit) f dt'TL{t')\L^,{t,t')\\ (11) 

J- Jo J -oo 

where Z/2i(t,t') is a matrix element of the operator L(t,t'); L2i{t,t') = {e2, L(t,t')ei). 

Since the probability amplitude for tunnelling is exponentially sensitive to the position of 
the dot, the maximum of the tunnel exchange interaction between an electron on the dot and 
an electron in one lead occurs when the tunnelling coupling to the other lead is negligible. 
This is why we will assume the following property of tunnelling amplitude Ai_ij(t) to be 
fulfilled: 

Ti(t)T^(t) = 0, Ti(t),T^(t)^0 (12) 

This assumption allows us to divide the time interval (0, T) into the intervals (0, r) + 
(r,r/2) + (T/2,r/2 + r) + (T/2 + r,T). We suppose that TL{t) ^ (but H = 0) only 
in the time interval (0,r) (and, analogously, Tji{t) = Tiit + T/2) in the time interval 
(T/2, T/2 + r)). Using this approximation together with the properties (fTUj) of the operator 
L[t,t'), we arrive at the following expression for the average current (Appendix 2): 

OD 

/ = |(l-e-T5^|(e2,L(T/2,r)L"e0p. (13) 

n=0 

Here L = L(T + t,t) and 

T = 2nu [ dtTl{t) (14) 
Jo 

is the tunnelling rate. Consequently, in order the calculate the average current it is necessary 
to investigate the properties of the evolution operator L. It follows from its definition that 

L = L{T + T, t) = L{T + r, T)L(T, T/2 + r)L(T/2 + r, T/2)L(T/2, r) (15) 
= e~^^+'^^^'"/^+*'^^*°L(T, T/2 + T-)e~^^~'^^^'"/'^"*'^^*''L(T /2, r) , 
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where $o = Jq dtJ{t) . From the symmetry properties of the operator i?(T/2 + r < t < T), 

= R, a2R* = -R, asRi-t) = -R{t)a3 
it follows that the operator U = L{T, T/2 + r) has the form 



U=\^'-'' I (16) 

In addition to this, L(T/2, r) = aiUai. As a result, the operator L can be expressed as 

L = (e-'^^^/^+^*«^^f/ai)' . (17) 

Proceeding with the analysis we (i) calculate the eigenvalues Aj and eigenvectors b, of 
the operator L of Eq. (fTTjl : Lbj = Ajbj, (ii) substitute the expansion = ajihj (where 
(a)~^ = (ej,bj)) into Eq. (fT^ and calculate the average current. The result is 



CK , ^, cosh r/2 + cos 2'(9 . ^ 

/=— sinhr/2 ^- / — , (18) 

T ' sinh^ r/2 + /€(1 + cos 2t9 cosh r/2) ^ ^ 



where 'd = if + ^o,k = 27^/(1 + 7^). Equation (fTHj) for the average current is our main 
result. 



CALCULATION OF THE CURRENT IN LIMIT OF STRONG AND WEAK 
EXCHANGE COUPLING BETWEEN THE DOT AND THE LEADS 

Although the result ()18j) for the tunnel current is both transparent and compact, it is in 
general a rather complicated problem to find the magnetic field dependence of the coefficient 
K, which depends on the probability amplitude 7 for flipping the spin of shuttled electrons. 
Three different time scales are involved in the spin dynamics of a shuttled electron. They 
correspond to three characteristic frequencies: (i) the frequency of spin rotation, determined 
by the tunnel exchange interaction with the magnetic leads; (ii) the frequency of spin rotation 
in the external magnetic field, and (iii) the frequency of shuttle vibrations. Different regimes 
occur depending on the relation between these time scales. Here we will consider two limiting 
cases, where a simple solution of the problem can be found. Those are the limits of weak 
•Jl{r) f^H and strong Jl{r) ^ fJ'H exchange interactions with the leads. 
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Weak exchange interaction 

In the limit Jl{r) f^H one may neglect the influence of the magnetic leads on the 
on-dot electron spin dynamics. In this case the matrix U given by Eq. pfi|) can easily be 
calculated and Eq. ([T^ reduces to 

2e sin2^9/2 tanhr/4 



T sin^ ^9/2 + tanh2r/4 

-T/2 



(19) 



where 'd = gn dtH is the rotation angle of the spin in the external field. 

Two different scales for the external magnetic field determine the magneto-transmittance 
in this limit. One scale is associated with the width of the resonant magnetic field dependence 
(see the denominator in Eq. ()19|) ). This scale is (restoring dimension) 

6H = T^, (20) 

where uj is the shuttle vibration frequency. The second scale, 

AH = ^. (21) 

91-^ 

comes from the periodic function sin^'(9/2 that enters Eq. p9|) . The magnetic- field depen- 
dence of the current is presented in Fig. 2a. Dips in the transmittance of width 6H appear 
periodically as the magnetic field is varied, the period being A. This amount to a giant 
magneto-transmittance effect. It is interesting to notice that by measuring the period of 
the variations in I{H) one can in principle determine the shuttle vibration frequency. This 
amounts to a dc method for spectroscopy of the nanomechanical vibrations. Equation ()2ip 
gives a simple relation between the vibration frequency and the period of the current varia- 
tions. The physical meaning of this relation is very simple: every time when u/Q = n + 1/2 
{Q is the spin precession frequency in a magnetic field) the shuttled electron is able to fully 
flip its spin to remove the "spin-blockade" of tunnelling between spin polarized leads having 
their magnetization in opposite directions. 



Strong Dot-Leads Exchange Interaction 



A strong magnetic coupling to the leads, Jmax ^ fJ'H, preserves the electron spin polariza- 
tion, preventing spin-flips of shuttled electrons due to an external magnetic field. However, 
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FIG. 2: Magnetic-field dependence of the transmittance of the device shown in Fig. 1 for the 

hmiting cases of a) weak and b) strong exchange couphng between dot and leads. The period AH 
and the width 5H of the "dips" are given by Eqs. (21) and (20) for case a) and 6H is given by 
Eq. (24) for the case b). 



if the magnetization of the two leads are in opposite directions, the exchange coupling to 
the leads have different sign. Therefore, the exchange couplings to the two leads tend to 
cancel out when the dot is in the middle of the junction. Hence the strong exchange interac- 
tion affecting a dot electron depends on time and periodically changes sign, being arbitrary 
small close to the time of sign reversal. In Fig. 3 the on-dot electronic energy levels for spins 
parallel and antiparallel to the lead magnetization are presented as a function of time. The 
effect of an external magnetic field is in the limit Jl{r) ^ fJ'H negligible almost everywhere, 
except in the vicinity of the level crossing. At this "time point" ,which we denote tiz, the 
external magnetic field removes the degeneracy and a gap is formed in the spectrum (dashed 
curve). The probability of electronic spin-flip in this case is determined by the probabil- 
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FIG. 3: On-dot energy levels for spin- up and spin-down electron states as a function of the position 
of the dot. Level crossing in the middle of the device is removed by an external magnetic field. 

ity of a Landau-Zener reflection from the gap formed by the magnetic field (in this case a 
Landau-Zener transition across the gap is a mechanism for backscattering of the electron, 
since this is the channel where the electronic spin is preserved). The matrix U can readily 
be expressed in terms of Landau-Zener scattering amplitudes. The amplitude and phase of 
electronic spin-flip is given hj ip = ipo + ^i,ipo is the Landau-Zener phase shift, 



$1 



T/2-T 



dtJit) 



and 7^ is the probability of the Landau-Zener "backward" scattering. 
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1 — exp 



" J'itLz) 



(22) 



(23) 



Schematical view of I{H) dependence is presented on a Fig. 2b. The width 5H of the 
minimum in I{H) dependence can be found directly from Eqs. (jl8|) . (j23|) 



5H 



(24) 



where Jo = mm{JL{R)(t)). 



CONCLUSION 



The analysis presented above demonstrates the possibility of a giant magneto- 
transmittance effect caused by shuttling of spin-polarized electrons between magnetic source- 
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and drain electrodes. The sensitivity of the shuttle current to an external magnetic field 
is determined, according to Eq. (pUj) . by the transparency of the tunnel barriers. By di- 
minishing the tunnelling transmittance one can increase the sensitivity of the device to an 
external magnetic field. The necessity to have a measurable current determines the limit of 
this sensitivity. In the low transparency limit, F <^ 1, the current through the device can 
be estimated as / ~ eTu. If one denotes the critical field that determines the sensitivity of 
the device by Her, one finds from Eq. (PUI) that Her — ^H. The critical field can now be 
expressed in terms of the current transmitted through the device as 

Her{Oe) ~ ~ ^(3 X 102)/(nA) , (25) 
e/i^ g 

where qq (= 2) is the gyromagnetic ratio for the free electrons. For / ~ 10^^ ^ 10^^ nA 
and qq/ g ^ 1/3 this gives a range H^r — 1 ^ 10 Oe. A further increase in sensitivity would 
follow if one could use a shuttle with several (A^) electronic levels involved in the tunnelling 
process. The critical magnetic field would then be inversely proportional to the number of 
levels, H,r{N) = H,r{N = 1)/N. 
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Appendix 1 

The Shrodinger equation results in equations for the coefficients c^'^j;^^ , c^^^ : 

^^ = -J{t)c^ - igf^H/2)c'f + n{t) ^ cf (t) , (26) 

t-^ = Jit)cf - {g^H/2)cf + Tnit)J2cf{t) , 

= e,cf + n{t)c^t), 
= e,cf + Tn{tX{t). 
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As it follows from the last two equations (together with the initial conditions) 

cf (t) = e-''^'5^p - I f dt'e''^^''-'^TL{t')c^t') , (27) 

J — oo 

cf{t) = -I f dt'd'^'^''''^TR{t')q{t'). 

J — oo 

Therefore, for the '^pC^it) one gets 



a J — oo fD 



/3 /3 

In wide-band approximation we suppose v{e) = const, therefore e*^'' = 2Txv5{t' — t) 
and 

Y^cf{t) = -t^uT^{t)c1. (28) 

Analogously, 

cf{t) = e-''"' - i7iuTL{t)c'f . (29) 

Substitute the expressions, Eqs. (|^Hj) . (j^ . to the first two equations (PBj) . one get the equa- 
tion Eq. (jH)) for the bi-vector c". 

Appendix 2 

Under our approximation we can change the integration limits in Eq. 

/ = i2nuy ^ £ dtTUt) f_ dt'Tl{t')\L,,it,t')\' 

pT/2+T pr 

= {27^vf- / dtTlit) / dt'Tl{t')\L,,{t,t')\\ (30) 
Jr/2 J-oo 

Beside this, in the time moments T/2 <t< T/2+t L(t, T/2) is a diagonal matrix. Therefore 
L2i(t,t') = i/22(t, T/2)I/2i(T/2, t'). As a consequence, the integral in the expression for the 
average current, Eq. dHHI), is factorized: 

I={27rur- / dtT'^{t)\L,,{t,T/2)\' / dt'Tl{t')\L,,{T/2,t')\' . (31) 

J-oo 

The first integral in Eq. (j31|) is easy to calculate. Having in mind that (T/2 < t < T/2 + r) 

ft 



|L22(t,T/2)|^ = exp 



27rz/ /" dtT^{t) 

Jt/2 
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one gets 

j-T/2+r 

2nu c/tT|(t)|L22(t,r/2)|2 = 1 -e"^, (32) 

Jt/2 

where quantity F is defined in Eq. (fT^ . 

The calculation of the second integral in Eq. ()31|) can be done in the same manner. One 
has the set of equahties, 

/T °° p~nT+T 

Tl{t)\L,,{T/2,t)\' = J2 / dtTl{t)\L,,{T/2,t)\' 

CO pj- 

J2 dtTl{t)\{e2,L{T/2,T)L{T,t-nT)e,)\\ (33) 



n=0 



For the quantity L{T,t — nT) = L{t + nT, t) one has 

L{T+nT, t) = L{T+nT, r+(n-l)T)L(r+(n-l)T, r+(n-2)T)...L(r, t) = U\t+T, t)L{t, t) . 

(34) 

Therefore, 

r dtTl{t)\L,,{T/2,t)\' = f2 r dtTlme,,LiT/2,T)L-L{T,t)e,)\\ (35) 

J-oo Jo 



Calculating the integral in the same manner, as in Eq. ()32j) . one gets the Eq. ()13j) for the 
average current. 
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